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ABSTRACT 

Time-series transit photometry from the Kepler space telescope has allowed for the discovery of 
thousands of exoplanets. We explore the potential of yet improved future missions such as PLATO 
2.0 in detecting solar system analogues. We use real-world solar data and end-to-end simulations to 
explore the stellar and instrumental noise properties. By injecting and retrieving planets, rings and 
moons of our own solar system, we show that the discovery of Venus- and Earth-analogues transiting 
G-dwarfs like our Sun is feasible at high S/N after collecting 6yrs of data, but Mars and Mercury 
will be difficult to detect due to stellar noise. In the best cases, Saturn’s rings and Jupiter’s moons 
will be detectable even in single transit observations. Through the high number (>1 bn) of observed 
stars by PLATO 2.0, it will become possible to detect thousands of single-transit events by cold gas 
giants, analogue to our Jupiter, Saturn, Uranus and Neptune. Our own solar system aside, we also 
show, through signal injection and retrieval, that PLATO .S.d-class photometry will allow for the 
secure detection of exomoons transiting quiet M-dwarfs. This is the first study analyzing in-depth the 
potential of future missions, and the ultimate limits of photometry, using realistic case examples. 
Suhjeet headings: planets and satellites: detection 


1. INTRODUGTION 

High precision, high duty time-series photometry from 
Kepler has contributed to numerous fundamentally new 
discoveries in the exoplanet field (e.g. Borucki et al. 
(2010); Burke et al. (2014)). The primary mission ended 
in 2014 after finding thousands of planets and candidates 
(Mullally et al. 2015), with the technical failure of two 
reaction wheels. The extended K2 mission (Howell et al. 

2014) continues to deliver data and new planets (Van- 
derburg et al. 2015; Foreman-Mackey et al. 2015). 
After this huge success, the next spacecraft photometry 
missions with improved technology are expected for 2017 
(Transiting Exoplanet Survey Satellite, TESS (Ricker et 
al. 2014)) and 2024 (Planetary Transits and Oscillations 
of stars, PLATO 2.0 (Rauer et al. 2014)). We asked 
the obvious question: What can we expect from these 
missions in the very best case? And ultimately: Assum¬ 
ing near-perfect photometric technology at some point 
in the future, what can we expect from photometry as 
such? Where are the fundamental limits? 

To begin, we will examine and define instrumental and 
stellar noise (section 2). Afterwards, we will discuss the 
Kepler, TESS and PLATO 2.0 mission designs and lim¬ 
itations (section 3). We will not focus on the high num¬ 
bers of discoveries, but on the very best cases with re¬ 
spect to their instrumental noise. In section 4, we will 
present the view of a distant observer at our planets 
transiting our Sun, assuming near-perfect photometry. 
Inversely, this is what we can expect from future space 
missions when it comes to finding solar system analogues. 

We will conclude with an outlook to the limits of pho¬ 
tometry. 

hippke@ifda.eu 
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2. EXAMINING THE NOISE 

Noise in data is often the limiting factor of data anal¬ 
ysis. Noise in exoplanet transit photometry is caused by 
instrumental imperfections (Nj) and stellar jitter {Ns). 
In the following, we will discuss both parts separately. 
The total noise N, assuming Gaussian distribution, is 
then calculated as: 

N=^Nf + Nl ( 1 ) 

2.1. Stellar noise characterization 

Stellar noise occurs with different characteristics on all 
time scales. First considerations for the Kepler mission 
by Batalha et al. (2002) were theoretical, due to the 
lack of precise data for other stars. In our sun, there is 
the 11-year solar activity cycle (Schwabe 1843; Usoskin 

2009) , a phenomenon we also find on different time scales 
in other stars (Garcia et al. 2010). The solar rota¬ 
tion of ^27 days (Bartels 1934; Beck 2000) introduces 
noise from spots, first noted by Galileo Galilei in 1612 
(Scheiner 2010). It has also been argued that the so¬ 
lar activity is modulated by planetary gravitational and 
electromagnetic forces acting on the sun, namely those 
by Mercury, Venus, Earth and Jupiter (Scafetta & Will- 
son 2013). 

In the following, we will ignore trends longer than a 
few days, and assume they can be removed using filters 
such as Savitzky & Golay (1964) used by Gilliland et 
al. (2011) to analyze Kepler noise, median filtering (e.g. 
Garter at al. (2012); Tal-Or et al. (2013)) or polyno¬ 
mial fitting (e.g. Santerne et al. (2014); Gautier et al. 

(2012)). Instead, we focus on the jitter on time scales 
of planetary transits, mostly 1-10 hours (Koch et al. 

2010) . This jitter originates mainly from stellar oscilla¬ 
tion modes, granulation at the surface of the star, and 
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Brightness (Kp) 

Fig. 1.— Kepler instrumental noise (straight line) depends on 
the brightness of the star. This plot shows only the brightest 
end, and only G-dwarfs. Our sun would exhibit a total noise 
(instrumental plus stellar) in between the dashed lines, which 
give the limits for the active (upper) and quiet (lower) sun. Few 
other stars are more quiet than our sun, and many are more active. 

rotational activity (Andersen & Korhonen 2015). 

Our sun’s noise varies by a factor of ^2 during the 11- 
year solar cycle, from 7.8ppm (2007.77, quiet period) in 
6.5hrs bins to 14.7ppm (2002.39, active period) (Gilliland 
et al. 2011; Frohlich et al. 1997). This is at the quiet 
side of G-type stars, of which the most quiet 1% have 
6ppm, with a total cut-off at 5ppm (Ghristiansen et al. 

2012; Basri et al. 2013). Although the noise mea¬ 
sures and results differ slightly between these authors, it 
can also be seen from Figure 1 that there are few stars 
more quiet than our sun. There is a strong dependence 
of stellar noise to temperature: Gooler stars are usu¬ 
ally more active, so that among M-dwarfs only very few 
are as quiet as our sun, and most are around ^50ppm. 
However, there might be a few extremely quiet (l-4ppm) 
G-dwarfs (Hall et al. 2007), theorized to exhibit a time 
of almost no spots, as was the case for our Sun dur¬ 
ing the Maunder Minimum (Maunder 1912; Zolotova 
& Ponyavin 2015). The detection of such a fortunate 
case, where mainly granulation (Ippm) contributes to 
stellar noise, would allow for extreme observations, given 
near-perfect technology. As we have not detected such 
a very-low noise star yet, we will instead concentrate on 
the known quiet end of G- and M-dwarfs. 

Noise in sun-like stars, on this timescale, originates 
mostly from asteroseismic oscillations. Solar-like oscilla¬ 
tions, mostly acoustic or pressure (p) modes, are reported 
with their highest amplitudes at frequencies between 
29.9/iHz (10 hours) and 3619^Hz (4 minutes) among Ke¬ 
pler stars (Huber et al. 2012). Regarding the solar 
noise properties on time scales of ~hours, it can clearly 
be seen from Figure 2 (top) that the jitter is not Gaus¬ 
sian. Glearly, there are trends, mostly from oscillations 
and spots, lasting a few hours, showing spots appearing 
and disappearing on the disc. This is the same timescale 
on which exoplanet transits occur, so that these trends 
cannot easily be filtered out. Time-correlated red (Brow¬ 
nian) noise becomes more and more Gaussian with stacks 
from different epochs, because the correlation of the total 
noise decreases (Barnes & Allan (1966) and references 
therein). This will be explained in more detail in sec¬ 



Fig. 2.— Data for the sun from VIRGO/DIARAD (Frohlich 
et al. 1997) shows trends on time scales of hours and longer 
(top). We inject synthetic transit shapes into these raw data 
(middle). To retrieve the transit, we fit a long-time sliding median 
while masking the times of transit (bottom). Instrumental noise 
(~0.4ppm) is smaller than the symbol size. 

tion 2.4. 

Another noise source is contamination of the starlight 
from background (or foreground) stars in the aperture, 
which can occur when the angular separation of the con- 
taminator is smaller than the spatial resolution of the 
instrument. For Keplers (low) spatial resolution, this is 
a problem in some cases, but not in general: “The overall 
increase of median and mean noise (...) are only 0.2 and 
0.1 ppm” (Gilliland et al. 2011). 

2.2. Stellar noise modeling 

Stars show brightness variations on different time 
scales, with different amplitudes, and characteristics. On 
the extreme side, there are for example RRab Lyrae with 
an amplitude of ~50% over 0.5d, but near-perfect repeti¬ 
tion in some cases (Smith 2004; Szabo et al. 2014). On 
the quiet side, there are stars like our Sun, or Kepler- 
197 (Rowe et al. 2014), which show very low (0.1%) 
long-time (months) variation, but stochastic trends on 
an hour to day timescale. Intermediate cases are most 
common, e.g. Kepler-264b (Rowe et al. 2014; Hippke 
2015) with strong short-time trends, or Kepler-96 (Marcy 
et al. 2014) with prominent long-time trends. We show 
these examples in Figure 3 and use them for the following 
discussion of suitable noise models. 

The most fortunate case is to have trends that are 
on a much longer timescale than the transit signature. 
Kepler-96 is an example for this behavior, where a “stel¬ 
lar noise model” using splines or polynomials will suc¬ 
cessfully remove the trends which are induced mostly by 
stellar rotation. The short transit times can be masked 
for the detrending, and the transit analysis can be done 
afterwards. A simultaneous fit would only complicate 
the process, and have no additional benefit. The limits 
of this method are reached in cases where the trends are 
less sinusoidal, as is the case for CoRoT-7b (Haywood et 
al. 2014; Barros et al. 2014). Then, Gaussian Processes 
(GPs) are suitable to treat these systematics. Their main 
advantage is to naturally handle correlations irrespective 
of their origin, by specifying high-level properties of the 
covariance (Evans et al. 2015). Gommonly used correla- 
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Fig. 3. — After removing instrumental trends, different charac¬ 
teristics in stellar noise are evident. Kepler-96 (top, black) shows 
strong long-time (weeks) variations, but is as quiet as the other 
examples on hours to days timescale. Kepler-197 (middle, red) is 
a star with low noise on all time scales, like our Sun. Kepler-264b 
(bottom, blue) is as quiet long-time, but shows more prominent 
short-time (hours, days) trends. All stars have been centered to 
exhibit a transiting planet at T=5d. 

tions are radial-velocity measurements (Haywood et al. 
2014), different wavelengths (Evans et al. 2015), or in¬ 
strumental systematics such as the drift in the roll-angle 
of Kepler’s K2 mission (Aigrain et al. 2015; Foreman- 
Mackey et al. 2015). 

When the noise timescale gets closer to the transit 
duration, other methods must be chosen. Most promi¬ 
nently, the wavelet-based formalism described in Carter 
& Winn (2009) is being used, e.g. by Huber et al. 
(2013b) in their analysis of spin-orbit misalignment in 
the multiplanet system Kepler-56, or by Barclay et al. 
(2015) in their confirmation of Kepler-91 being a “giant 
planet orbiting a giant star”. Using wavelets, the noise 
is modeled by the sum of two stationary (Gaussian) pro¬ 
cesses; one is uncorrelated in time, and the other has a 
spectral power density as 1/P ■ 

In any such model, parameters must be estimated and 
validated individually for each star, using Monte-Carlo 
simulations. This is also true for competing methods, 
such as the “time-averaging” method (Pont et al. 2006), 
or the “residual-permutation” method (e.g., Jenkins et 
al. (2002)). Consequently, stellar noise modeling is 
not applied to the large number (thousands) of planets 
found with automatic pipelines. For interesting individ¬ 
ual cases, however, stellar noise modeling can and should 
be used. It is important to note that these cases have, 
beforehand, been detected without modeling, so that any 
detection must pass this initial threshold^. In the course 
of this paper, we choose the same approach: We use the 
“standard” pipeline without noise modeling, but intro¬ 
duce it after the detection of the Earth, which can be 
assumed to be of distinguished interest for any distant 
observer. Noise modeling reaches its current limits for 
signals so small that they become indistinguishable from 
permanent noise features; for our Sun, this would be the 
case for signals on ppm-level as for Earth’s moon. There- 

^ Further advantages in computer algorithms might introduce 
methods that can self-adjust to individual cases and validate the 
result. 


fore, we limit our example to Earth itself. 

2.3. Instrumental noise 

Even the most perfect instrument will produce some 
noise. Fundamentally, this originates from the fact that 
photons (starlight) and electrons (detector) are quan¬ 
tized (Einstein 1905), so that only a finite number can 
be counted in a given time. This phenomenon is the shot 
noise (Schottky 1918), which is correlated mostly to the 
brightness of the target. In addition, noise occurs from 
the readout of the CCD, when the small signal gets am¬ 
plified. On a timescale of 6.5hrs (13 bins of ^30min), the 
Kepler instrumental noise for a Kp=12 star from Pois¬ 
son (shot) and readout is 16.8ppm. Two other instru¬ 
mental noise sources have been quantified for the Kepler 
spacecraft: Intrinsic detector variations (lO.Sppm), and 
a quarter-dependent term (7.8ppm) (Gilliland et al. 
2011). The total instrumental noise is then 20.4ppm. 
For reference, the design of the spacecraft expected in¬ 
strumental noise of 17.0ppm. 

2.4. Signal-to-noise definition and noise 
characterization 

It is a common practice in exoplanet science (Jenkins 
et al. 2002; Rowe et al. 2014) to define the signal-to- 
noise ratio as the depth of the transit model, compared 
to the out-of-transit noise: 

S/N=V^— (2) 

aoT 

with Nt as the number of transit observations, T^ep the 
transit depth and gqt the standard deviation of out-of¬ 
transit observations. This measure overestimates S/N 
for large planet-to-star radius ratios, and large impact 
parameters, but these configurations will be neglected 
in the present work. As pointed out by Fressin et al. 
(2013), the detection of KOIs becomes unreliable for a 
S/N < 10. This assumes Gaussian noise - the noise dis¬ 
tribution of our Sun, even during the quiet year 2007, 
is distinctly not Gaussian. It has strong time-correlated 
(red) noise features and a “long tail”. Two statistical 
methods should be employed to quantify this. First of 
all, we recommend a skewness and Kurtosis tests for nor¬ 
mality, e.g. by Shapiro & Wilk (1965). For our sun, as 
expected, the difference to Gaussian noise is significant, 
at the 0.1% level. After establishing that the noise is not 
Gaussian, we propose to measure the “redness” per bin 
length L, as described by Steves et al. (2010): 

fTbinned = (^unbinned (3) 

For completely uncorrelated noise, /? = —1/2, while for 
Gaussian noise /3 = 0. Instead of giving ,0-values for 
each bin length, we chose to calculate the achieved per¬ 
centage (where Gaussian noise is 100%), and show the 
result in Figure 4. This gives the recoverable S/N of 
transits, when compared to Gaussian noise. As shown 
in the figure, the S/N in solar red noise is significantly 
lower. This issue originates from the simple fact that 
subsequent data points are time-correlated, and are thus 
more likely to have (nearly) the same error (positive or 
negative) when binned together on the time axis. In 
other words: Red noise data doesn’t bin as well as white 
noise data. The penalty of this noise characteristic is 
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Fig. 4.— S/N ratio for transits that can be achieved with real red-noise solar data, when compared to theoretical pure white Gaussian 
noise. Left: Observed (data points, straight line) and simulated Gaussian noise (dashed line). Right: O-C for the noise. The structural 
break at 7hrs is evident. For Earth-long transits, the achieved S/N is 49%, and higher for longer transit durations. 


also time-correlated, because, over time, it diminishes. 
For an Earth-analogue transit (13.2hrs), only 49% of the 
expected S/N can be achieved; 60% for a Jupiter (33hrs) 
and 71% for a Neptune (76hrs) (Figure 4). When stack¬ 
ing different years (of similar noise), the penalty is zero, 
as there is no time-correlation any more. The original 
intra-transit penalty cannot be recovered, of course, but 
stacking different epochs brings the expected (like Gaus¬ 
sian) stacking bonus. 

As shown in Figure 4, the expected power law when 
binning data (for mainly Gaussian noise) begins at 
^7hrs. For shorter times, the best-fit is given by a lin¬ 
ear correlation (i?^ = 0.99; best-fit power law gives lower 
E? = 0.87). The location of the structural break can be 
calculated using the test by Chow (1960), asking the 
question whether the coefficients in two regressions for 
different data sets are equal. The result is a clear struc¬ 
tural break at the 0.1% level, with a best-fit location at 
7±lhrs. To sum up, our sun produces Gaussian noise on 
long (l^ld) timescales, but suffers red noise punishment 
on shorter timescales. 

3. METHOD: PREPARING THE TELESGOPES 

In this section, we will discuss the instrumental per¬ 
formance of the future TESS and PLATO 2.0 missions. 
We will define the properties for two “telescopes” to use: 
The best-case PLATO 2.0 performance, and a theoret¬ 
ical near-perfect instrument, dubbed PERFECT. These 
instruments will then be used in section 4 to inject and 
retrieve transits. 

3.1. Observations with TESS 

The longest observing intervals for TESS (Ricker et 
al. 2014) during its 2-year mission will be one year (for 
parts of ths sky near the ecliptic pole). It will focus on the 
brightest {Ic = 4—13) stars, suitable for follow-up obser¬ 
vations with the James Webb Spaee Telescope (Gardner 
et al. 2006) (the areas of the sky with one year coverage 
from TESS will also fall into the continuous viewing zone 
for JWST) . As such, its mission design is not targeted 
at the very best photometric performance, although it is 
expected to be on par with Kepler (Ricker et al. (2014), 
their Figure 8). In this paper, we strive for best possible 


photometry, and are interested in solar system analogues, 
so that we will neglect TESS, as it does “not address the 
science case of characterizing rocky planets at intermedi¬ 
ate orbital distances (a>0.3au, including the HZ) around 
solar-like stars, which remains unique for PLATO 2.0.” 
(Rauer et al. 2014). 

3.2. Observations with PLATO 2.0 

Instrumental noise from PLATO 2.0 is expected to be 
as low as lOppm on 30min timescale for bright (V=9) 
stars (Zima et al. 2006), or 8ppm in Ihr integrations 
(Rauer et al. 2014), when observed by many (up to 
36) cameras simultaneously, thus producing ^3ppm of 
instrumental noise on the same 6.5hrs timescale. While 
Kepler observed ~30 bright (V<8) stars, PLATO 2.0 
will collect data for ~3,000 - giving a good chance for a 
useful share of quiet stars among them. Observing our 
sun through PLATO 2.0 would then give a stellar noise 
fraction of 74% (quiet sun) to 84% (active sun), making 
such an observation strongly limited by intrinsic stellar 
noise. For smaller, but generally more active M-dwarfs 
(with noise levels around SOppm (Basri et al. 2013)), 
this ratio can reach 95% to 99%, making the technology 
near-perfect and observations only limited by intrinsic 
stellar noise. This shows the ultimate limits of photom¬ 
etry as such, and the importance of understanding (and 
modeling) stellar noise. PLATO 2.0 has an expected 
lifetime of 6 years, giving a useful duration for multiple 
transits of Earth and Mars analogues. 

Eor simulating realistic data as expected from PLATO 
2.0, an end-to-end simulator^ evolved over the last 
decade (Zima et al. 2006; Marcos-Arenal et al. 2014). 
It takes into account effects down to the sub-pixel ma¬ 
trix, satellite orientation jitter, PSF convolution and all 
CCD-related noise sources. Using this simulator, we have 
created a flat lightcurve for a bright (V=9) star, resulting 
in Gaussian noise at ~10ppm in 30min bins. These data 
are the instrumental basis of our PLATO 2.0 simulations 
in section 4. 

3.3. The PERFECT telescope 

^ https://fys.kuleuven.be/ster/Software/PlatoSimulator 
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Fig. 5.— Transit light curves for solar system inner (left) and outer planets (right, with Earth for comparison). Impact factor has been 
set to 0, so that the transit duration is maximized. These artificial curves are used for the injections. 


PLATO 2.0 is expected to launch in 2024 and end its 
nominal operation in 2030 For the following decade 
(the 2030s), one might imagine a successor mission, with 
even better instruments. We should assume that many 
noise sources can be reduced (or eliminated completely) 
using the knowledge and expertise gained throughout 
previous missions, such as pointing jitter and thermal 
variations. Sensor sensitivity, for instance, has improved 
over the decades and approaches >80% quantum effi¬ 
ciency today, over the wavelength range 0.6-2.5/im (e.g., 
McGurk et al. (2014). 

For those stars where the photon flux (shot noise) is 
not the limiting factor, instrumental noise will be dom¬ 
inated by these other sources. As we have seen, total 
instrumental noise from Kepler is at ^20ppm, whereas 
we expect ~10ppm from PLATO 2.0. It will be interest¬ 
ing to check the perfect telescope with zero instrumen¬ 
tal noise. As explained in section 2.3, this is in princi¬ 
ple unphysical, but a noise floor down to 0.4ppm as for 
VIRGO/DIARAD might be achievable. We will keep 
(adopt) this 0.4ppm of instrumental noise for the virtual 
PERFECT telescope. In section 4.2, we will see that the 
difference between PLATO 2.0 and PERFECT is negli¬ 
gible for the standard quiet G-dwarf, as stellar noise is 
the dominating noise source. For a very quiet M-dwarf, 
however, this can make an improvement of up to 50% in 
S/N. 

3.4. Target stars and planetary systems 

As host stars, we choose a very quiet G2-dwarf like 
our sun - in fact, we will simply use our sun’s data from 
VIRCO/DIARAD (Frohlich et al. 1997; Appourchaux 
et al. 1997). Its bandpass is comparable to the pho¬ 
tometry space telescopes, and its instrumental noise is 
<0.4ppt and can thus be neglected. The performance 
is better than Kepler or PLATO 2.0 due to the high 
flux the instrument can receive from our nearby Sun, 
when compared to stars many parsecs away. We have 
interpolated the data to 30min bins (analog to Kepler 
long cadence data) and added aforementioned ^lOppm 
instrumental noise to these data in order to simulate a 

^ http://wi™2.Warwick.ac.uk/fac/sci/physics/research/ 

astro/plato-science/pre-launch/ 
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Fig. 6.— Noise characterization of KIC7842386, the most quiet 
M-star out of 1015 (for which CDDP data is available) from Kepler^ 
with stellar noise CDDP=6.7ppm on 6.5hrs timescale. The dashed 
line represents (theoretical) pure Gaussian noise on a straight line 
in this log-log plot. Measured noise in bins (dots with error bars) 
is consistent with Gaussian noise. 

best-case observation by PLATO 2.0. 

The other reference star is a O.Si?© Ml-dwarf on the 
very quiet end of the distribution, exhibiting 7ppm of 
Gaussian (not time correlated) noise. From 1015 char¬ 
acterized Kepler M-stars, two exhibit stellar (GDDP, 
Christiansen et al. (2012)) noise <7ppm. We have 
taken the most quiet one, KIC7842386 (6.7ppm on 6.5hrs 
timescale) and analyzed its noise characteristics. As can 
be seen in Figure 6, the noise is Gaussian within the er¬ 
rors. While the number of quiet M-stars might be low, 
we take the noise level of this fortunate example for our 
injections. 

Due to their smaller radii, Ml-dwarfs are particularly 
suitable to observe transits. Their absolute luminosity is 
smaller at ~ 3.5% of that of G-dwarfs, so that they need 
to be closer to the observer by a factor of Lq /Lj to have 
the same apparent brightness for the same PLATO 2.0 
instrumental noise. 

Using these virtual instruments, we will observe dif¬ 
ferent bodies transiting the G2-dwarf and the Ml-dwarf. 
We inject planet transits following the standard Man- 
del & Agol (2002) model, including quadratic limb- 
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darkening with stellar metallicity (Claret & Bloemen 
2011), as implemented by PyAstronomy^. For the G2- 
dwarf, we observe the planets of our own solar system, 
exhibiting a wide range of transit depths and -durations 
(Figure 5). We set the impact parameter to zero, in or¬ 
der to maximize the transit duration. Special focus is 
on our Earth, including Earth’s moon (section 4.2), and 
the Jupiter system (section 4.4). Finally, we will explore 
Saturn, including the flux gain caused by the forward¬ 
scattering of its rings (section 4.5), as well as Uranus 
and Neptune (section 4.6). To sum up, this section asks 
the question what we would see of our own solar sys¬ 
tem if we were placed somewhere else in the galaxy, with 
near-perfect photometric equipment. 

For the Ml-dwarf, we selected to test a 2 .Oi ?0 Super- 
Earth with Ganymede-sized (O. 4 i? 0 ) moon (section 4.7). 

For all cases, we will use SOmin integrations as the 
shortest bin available. This is comparable to Kepler's 
long-cadence (LC) bins, and guarantees sufficient sam¬ 
pling for all transits in this work. Only for Saturn’s rings 
we find that a finer time resolution would be (marginally) 
beneficial. 

3.5. Injection and retrieval of transits 

We use the raw data from VIRGO/DIARAD (Frdhlich 
et al. 1997) and inject the synthetic transit shapes as 
described in section 3.4 and shown in Figure 5. An ex¬ 
emplary step-by-step injection and retrieval is shown in 
Figure 2, here neglecting instrumental noise for clarity. 
To retrieve the transits, we apply a sliding median with 
boxcar length of ~2 days, while masking the data points 
affected by the transit. For the longer transit durations 
of the outer planets, the boxcar length was increased ac¬ 
cordingly, leaving slightly larger residuals. This method 
has the advantage that we can use solar data with negli¬ 
gible (<0.4ppm) instrumental noise, but fully preserving 
stellar noise. To these data, we add instrumental noise 
as described in section 3.4. This method allows us to 
control and modify noise sources separately. 

3.6. Transit probability 

The transit probability of any planet is low, so that 
many stars need to be observed, in order to collect a 
useful sample (Kepler observed >100,000 stars). For cir¬ 
cular orbits, the transit probability can be calculated as 
(Borucki & Summers 1984): 

Ptr = — (4) 

a 

where R^, is the radius of the host star and a is the semi¬ 
major axis of the planet’s orbit. This gives transit prob¬ 
abilities e.g. for Earth as 0.47%, and Jupiter 0.041%: 
In order to potentially detect one Earth-analogue, one 
must on average survey 213 G-dwarfs that host an Earth 
like planet, and 2,439 that host a Jupiter. Assuming 
?70 = 0.1 (Burke et al. 2015), one must survey 2130 
G-dwarfs to detect an Earth-analogue. Therefore, ob¬ 
serving a large number of (dimmer and dimmer) stars 
is required at high sensitivity. The instrumental per¬ 
formance of PLATO 2.0 is expected to be sufficient for 
the detection of ~ 150 Earth-sized planets on 365d-orbits 
(Rauer et al. (2014), their Figure 5.5). 

https : //github. com/sczesla/PyAstronomy 





Fig. 7.— Transits of Mercury, Venus, and Mars with six years of 
PLATO 2.0 coverage and a quiet sun. Due to the decreasing transit 
frequency, the error bars get larger in this order. Transits of Venus 
are clearly detected, but Mercury and Mars are not recoverable. 
For better comparison, all axes are identical and the same as for 
the Earth transits in Figure 8. 

4. RESULTS: OBSERVING TRANSITS WITH PLATO 2.0 
AND PERFECT 

4.1. Mercury, Venus and Mars 

Mercury is the smallest (O.38i?0) planet in the solar 
system with a transit depth of only 13.1ppm. On the 
other hand, its short period of 87.97 days allows for the 
greatest number of observed transits for a given observa¬ 
tion time. When collecting six years of PLATO 2.0 data. 
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and no data loss occurs during transit, a total of 25 tran¬ 
sits can be recorded. We expect a nominal S/N=7.0 from 
this stack, but as can be seen in Figure 7, this is not suffi¬ 
cient for a detection due to the red noise characteristics. 
We conclude that Mercury analogues will likely not be 
found with photometry around normal G-dwarfs, when 
neglecting noise-modelling. 

The next planet, Venus, is somewhat more interest¬ 
ing, as it is roughly (O.Obi?®, 80.6ppm) Earth-size, and 
transits more frequently due to its smaller orbit. On 
the other hand, its transit duration is also shorter than 
Earth’s (0.46d vs. 0.55d), giving roughly similar detec¬ 
tion S/N (Venus: 23.8; Earth: 28.6) for both when using 
six years of PLATO 2.0 observations. Venus orbits in¬ 
side of the inner edge of the habitable zone (Kopparapu 
et al. 2013), but if orbiting a slightly less luminous stars, 
it might be habitable. 

Mars detectability, then, suffers from smaller size 
(O.53i?0, 25.4ppm) and longer orbit (only 3 transits 
in 6 years). Its transit duration (0.67d) helps only 
marginally with the transit detection, giving an insuf¬ 
ficient S/N=7.2, comparable to Mercury. 

4.2. Earth and moon 

As explained in section 2.1 and shown in Figure 2, red 
noise from stellar variation prevents the secure detection 
of an Earth-sized planet around a quiet G2-dwarf in a 
single observation, no matter how good the photometer 
is. Also, single transit-like events can occur from sudden 
instrumental sensitivity drops, as described for the false¬ 
positive moon of Kepler-90g (Kipping et al. 2015). The 
situation improves when stacking a few transits. For an 
expected mission duration of six years for PLATO 2.0, 
six transits could be observed of an Earth-analogue in 
the best case. As can be seen in Figure 8, the detection 
is then possible at high confidence for both the active 
(S/N=15.5) and quiet (S/N=28.6) sun. 

This finding is in agreement with the results of Aigrain 
& Irwin (2004), who found that a single Earth tran¬ 
sit cannot be detected, but a stack of six transits finds 
the planet at “high” S/N. While this result was achieved 
from synthetic data, Carpano & Fridlund (2008) have 
also used VIRGO/DIARAD data and concluded that 
planets smaller 2R^ cannot be detected in a single ob¬ 
servation. This shows that our method is valid and con¬ 
sistent with previous findings. 

The situation is equally promising for radial-velocity 
measurements (RV), assuming favorable conditions (sin¬ 
gle planet edge-on in circular orbit). As discussed by 
Lagrange et al. (2010) and Meunier et al. (2010), 
an Earth analogue can be detected with 1% false alarm 
probability after collecting > 200 epochs spread over 4 
years. 

Eor the detection of Earth’s moon, however, prospects 
are much worse. The tiny dip (6.6ppm) caused by Luna 
(i?(j = O.27i?0) is invisible in a single transit. The tran¬ 
sit timing variation (TTV) amplitude is 112 seconds, 
and the transit duration variation (TDV) only 14 sec¬ 
onds (Kipping et al. 2009). Errors from Kepler are on 
the order of 400 to 800 seconds (S/N^O.3), making a de¬ 
tection of this configuration impossible. As TDVs, and 
their detection sensitivity, are a function of the impact 
parameter, there might be a few cases that are more for¬ 
tunate. As explained by Kipping (2009), for very high 


impact parameters b > 1, i.e. grazing transits, TDVs 
are several times higher, up to 350s for realistic cases 
(see their Figure 2). Such candidates are hard to find, 
as the V-shaped transit curves are usually rejected by 
the automatic routines as eclipsing binaries. This pos¬ 
sibility however, together with shorter integration times 
to improve timing sensitivity, might make the detection 
of exomoons, based on timings, possible with PLATO 
2.0. It seems crucial to have integration times as low 
as technically possible (ideally less than a minute), and 
improvements in algorithms are required to not reject all 
V-shaped transits. 

Without noise modeling, finding an Earth’s moon ana¬ 
logue is also difficult using the orbital sampling effect 
(OSE) as first described by Heller (2014). In short, 
when adding up many randomly sampled observations, a 
photometric flux loss appears in the phase-folded transit 
light curve, reflecting the moon’s blocking of light. The 
effect depends mostly on the moon’s radius and planetary 
distance. The OSE can be used to detect a significant 
flux loss before and after the actual transit (if present), 
which might be indicative of an exomoon in transit. The 
basic idea is that at any given transit the moon(s) must 
be somewhere: They might transit before the planet, af¬ 
ter the planet, or not at all - depending on the orbit 
configuration. But by stacking many such transits, one 
gets, on average, a flux loss before and a flux loss after 
the exoplanet transit. While the OSE has shown to be 
useful when stacking many transits (Hippke 2015), the 
sheer number of ~180 transits required in this case (for 
a 2cr detection) renders the method useless for Earth and 
moon (Figure 9). 

4.3. Earth transit with noise modeling 

To explore the benefits of a tailor-made noise model¬ 
ing, we have implemented a standard wavelet-based filter 
as presented by Garter & Winn (2009) and explained in 
section 2.2. We simulate our Sun’s noise as a stationary 
process plus a time-correlated process of spectral power 
density 1//'’'. To estimate the best global parameters, 
we used Monte-Carlo runs to minimize the post-model 
squared residals, but also restricted the parameters to 
not fit trends on times < 2hrs in order to avoid over¬ 
fitting. The danger of such noise models is that with 
enough and sufficiently small parameters, one can fit out 
anything perfectly, be it a small transit (e.g. a moon), 
or a spot. For solar data, we find that the noise model 
is robust to filter out trends on timescales longer than 
a few hours, and more than a few 10s of ppm, so that 
an Earth-like transit signal benefits from noise modeling. 
Blindly recovering Earth’s moon, however, doesn’t work 
even with such a noise model, as the moon’s transit sig¬ 
nal (6.3ppm) is either still buried in noise, or has been 
removed altogether by too aggressive over-fitting. 

We show the application of our noise model in the bot¬ 
tom left panel of Figure 8. The uppermost (red) symbols 
show the raw data without an injection, and the simu¬ 
lated (modeled) wavelet-filter as the black line behind the 
datapoints. After injection of the transit signal, we have 
then performed a simultaneous fit to the transit curve 
and noise-model (blue symbols). The residuals (green) 
are shown in the bottom; these can be compared the the 
raw data (red). We have used the same solar data as in 
Figure 2, so that one can visually compare the benefit of 
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Fig. 8. — Earth observations after collecting data of six transits. Upper left: Kepler observing the Kp=12 quiet sun (total noise: 
26ppm). Middle left: The same observation for PLATO 2.0 (16ppm). Upper right: PLATO 2.0 with the loud sun (29ppm). Middle right: 
PERFECT observing the quiet sun (15ppm). The bottom panels show the result including our noise model, as explained in section 4.3. 
The blue symbols on the bottom left, and the panel on the bottom right, are created using a simultaneous fit of noise model and transit 
curve. The value for the quiet sun differs slightly from section 2.3 at this specific transit time. For better comparison, all axes are the same 
as in Figure 7 for Mercury, Venus and Mars. 
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Time from mid-transit (days) 

Fig. 9.— Orbital sampling effect (dashed line) for Earth’s moon 
using PLATO 2.0 and the six quiet years of the solar cycle (2005- 
2010). Luna cannot be recovered due to (mainly stellar) noise. To 
achieve a 2(7 detection, ~180 transits would be required. 

the noise model. 

In the bottom right panel of Figure 8, then, we zoom 
into a single transit, offering ~ 2x improved S/N. This 
shows the benefit of a noise model: It removes most of 
the time-correlated noise, and increases S/N at the same 
time. However, it is also suspect to removing smaller sig¬ 
nals, such as moons. Also, it requires considerable effort 
to implement a model, validate its parameters, and apply 
it to the data. In the future, we must hope to achieve 
further theoretical advantages so that such models can 
be created automatically during planet-searches. 

4.4. Jupiter and moons 

A single Jupiter transit produces a deep (1.125% = 
ll,250ppm), long (33hrs) transit dip, resulting in a 
highly significant (S/N=4,360) detection. Due to its long 
period of 11.86 years, we can not expect to observe a sec¬ 
ond transit with PLATO 2.0 during the spacecraft’s ex¬ 
pected 6-year lifetime. Assuming that the single transit 
itself is spotted in the data, we can search for accom¬ 
panying moons. In the case of Jupiter, all of its larger 
moons are almost sky-coplanar and would thus also be 
transiting. In this example, we examine the case of all 
moons being at maximum separation, as this case is usu¬ 
ally expected to yield the highest detection probability 
(Hippke 2015). For our Jupiter, this is not the case: due 
to the large (25i?p) separation of Callisto, this moon fin¬ 
ishes its transit (just) before lo goes into transit; there 
is no cumulative dip of all moons transiting simultane¬ 
ously (when at maximum separation). We neglect all 
other small moons, as they contribute less than 1% of 
additional transit depth. 

As can be seen in Figure 10, the whole transit ensemble 
can be detected with marginal confidence at S/N=8.6. 
Individual moons are not discernible, and it cannot be 
decided whether a single (large) moon, or a multiple 
moon configuration is observed. Such an observation 
would constitute a strong moon candidate worth follow¬ 
up observations, but likely no clear detection. 

4.5. Saturn with rings 

In this section, we examine our Saturn transiting our 
quiet Sun observed through PLATO 2.0. A simplified 


test for this configuration has been performed by Tusnski 
& Valio (2013), using synthetic Kepler data and a dark 
Saturn-model. In the following, we use the real solar data 
and the model from Barnes & Fortney (2004), which 
includes diffractive forward-scattering (Dyudina et al. 
2005). For our Saturn, data are available from the 1989 
occultation of 28 Sgr by Saturn (French & Nicholson 
2000), to adjust the model to the rings’ complex nature 
(see Figure 8 in Barnes & Fortney (2004)). The main 
difference to the pure black transit is the flux gain due 
to scattered light. We inject these data, together with 
the standard transit shape, and show the result in our 
Figure 11. We neglect any moons, as they have been 
treated in the previous section. As can be seen in the 
graph, the rings are clearly recovered, with the flux gain 
prominently seen. As can be seen in the figure, a time 
resolution of at least 30min is required, to avoid smear¬ 
ing of the light curve features. The out-of-transit fea¬ 
tures alone account for a S/N=21.5, fully sufficient to 
claim a detection without modeling. Full modeling is 
however encouraged, as ^50% of the flux delta occurs 
during transit, and most of it during ingress/egress. As 
pointed out by Zuluaga et al. (2015), rings cause an 
increase in transit depth that “may lead to misclassifi- 
cation of ringed planetary candidates as false-positives 
and/or the underestimation of planetary density”. By 
comparing results from astrodensity profiling (Kipping 
2014) to those from astroseismology (Huber et al. 2013), 
such anomaly low density planets could be detected. 

Saturn itself produces a highly significant S/N=3,575 
in a single transit. 

4.6. Uranus and Neptune 

Uranus and Neptune produce deep (l,382ppm, 
l,354ppm), long (2.45d, 3.07d) transits. The S/N for 
Neptune (783 vs. 714) is slightly higher, due to its longer 
transit duration. As can be seen in Figure 12, the tran¬ 
sit is visually compelling, but stellar noise may mimic 
moons where there are none. For a test, we have in¬ 
jected Neptune’s largest moon, Triton (O.21i?0, 4ppm), 
which is unrecoverable in the noise. Furthermore, Tri¬ 
ton is usually not seen in transit, due to its inclination 
of 129.6° (Agnor & Hamilton 2006) with respect to the 
Laplacian plane of the solar system. In such cases, even 
the largest moons would remain unnoticed. 

4.7. Super-Earth and moon transiting the M-dwarf 

We have now completed the tour through our solar 
system, and have seen that detecting planets and rings 
is possible for many scenarios. Exomoons, however, are 
more difficult to observe, making a detection unlikely 
for solar system analogues, as explained in sections 4.2 
and 4.4. We will thus try a more promising configura¬ 
tion: A very quiet M-dwarf (0.57?©, stellar noise: 7ppm 
CDDP), assumed to have a transiting planet of Super- 
Earth size (2.07?®) in a P=23.97d orbit, with an accom¬ 
panying Ganymede-sized (0.47?®) moon. This system is 
inspired by Kepler-236c and shares its stellar and plan¬ 
etary size, and their separation. The stellar noise prop¬ 
erties are as described in section 3.4, resembling the M- 
dwarf KIC7842386 with 6.7ppm Gaussian noise on 6.5hrs 
timescale. As can be seen in Figure 13, a single transit 
is clearly detected for the planet, but not for the moon. 
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Fig. 10. — Left: Transit configuration of Jupiter system, with all moons at maximum separation on the ingress side. When Ganymede is 
at mid-transit, Jupiter has not started its transit yet, and Europa has already completed its own. Right: Flux for this configuration (line). 
Dots are real data LC bins (SOmin), data points with error bars binned in 6.5hrs. Note how the Callisto egress occurs before planetary 
ingress, and before lo’s ingress. At no time, all four moons are transiting simultaneously, when at maximum separation. 






Fig. 11. — Upper left: Saturn transit without ring (black line) and theoretical ring solo transit (grey line). Upper right shows zoom. 
Higher than nominal flux is caused by diffractive forward-scattering of the rings. Lower left: Saturn including ring (black line) and data 
for PLATO 2.0 during quiet sun (errors as dot size). Lower right: Zoom into egress, where the ring is clearly detected. Data points are 
30min LC bins. 
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Fig. 12.— Uranus and Neptune transits with data as 30min LC 
bins. Neptune’s transit is slightly shallower but longer. Just before 
Neptune’s transit, a spot-induced 105ppm flux-drop occurs, which 
could be mistaken for a exomoon. 

With PLATO 2.0 photometry and 6 years of data (100 
orbits), the moon dip is clearly retrieved through its or¬ 
bital sampling effect (middle plot). The gain from PER¬ 
FECT (bottom plot) is in this case ~50%, thanks to the 
high share of instrumental noise (7ppm stellar, lOppm 
PLATO 2.0). Interestingly, PLATO 2.0 photometry al¬ 
lows for the detection of the “cap” in the OSE. This 
feature represents the relative flux gain before and af¬ 
ter planetary transit, when a large-orbit moon goes into 
stacked ingress (egress) before (after) the planet tran¬ 
sits. Detecting this feature gives more of a “shape” to 
the OSE, in contrast to a pure instrumental decrease 
in luminosity. To give an impression of the unfolded 
photometry, we show a riverplot (Carter at al. 2012; 
Nesvorny et al. 2013) in Figure 14. With 100 transits 
(each row is one period) of PERFECT data, the moon 
transit signature is visually evident even without folding. 
With PLATO 2.0 data, it is also visible when stretching 
the gray- scale. Such plots can be used for a sanity check 
of moons with PLATO 2.0 data. 

4.8. Oblate planets 

Planets are never perfect spheres, but show oblateness 
to varying degrees. Measuring the oblateness of exoplan¬ 
ets will help our understanding of planetary formation, 
rotation, and internal structure (Carter & Winn 2010). 
The transit light curve of an oblate planet is different 
from that of a spherical, with the main delta during 
ingress and egress. Zhu et al. (2014) have performed a 
search in Kepler data and find the hot-Jupiter HAT-P-7b 
to be a good candidate. Their analysis shows that Ke¬ 
pler-level photometry is sufficient to detect “Saturn-like 
oblateness (/ = 0.1) for giant planets {Rp/R^ = 0.1) 
around relatively bright (12 mag) stars”. 

To analyze the prospects of PLATO 2.0 in this re¬ 
gard, we have injected a Jupiter-sized planet including 
Saturn-like (/ = 0.1) oblateness with a 87-day (Mer¬ 
cury) period into 6yrs of quiet solar data. The calculated 
oblateness delta is the one from Zhu et al. (2014) and 
has kindly been provided by the author (Wei Zhu 2015, 
priv. comm.). We chose a configuration that generates 
a large oblateness signal, using a projected planet spin 
obliquity of 45°. As can be seen in Figure 15, the signal 
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Fig. 13.— 2.0/20 planet and Ganymede-sized (OAR^) moon or¬ 
biting a O.5i?0 M-dwarf with P=23.97d, resembling Kepler-236c, 
but with low noise properties (stellar noise CDPP=7ppm). Top: 
Single transit with PERFECT giving a clear dip for the planet 
(S/N=22.8), but low S/N=3.9 for the moon at maximum separa¬ 
tion. Middle: After stacking 6yrs (100 orbits) with PLATO 2.0, 
the orbital sampling effect clearly recovers the moon. The gain 
from PERFECT is in this case ~50% (bottom). 

is clearly retrieved, at S/N=16. As for planetary tran¬ 
sits, the signal-to-noise ratio will be higher for shorter 
periods, larger oblateness, and larger planet-to-star ra¬ 
tios. On the other hand, it will be harder to detect plan¬ 
ets with smaller oblateness, longer periods, and smaller 
radii (for a constant projected planet spin obliquity). Us¬ 
ing the large sample expected from PLATO 2.0, we can 
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Fig. 14.— River plot of moon injection, compare Figure 13, using 
the PERFECT data. Dark colors are lower flux levels and indicate 
the transits. The horizontal axis of each plot represents one transit 
event, with the vertical dark stripe being caused by the planetary 
transit. Each of the 100 rows is one period. Top: Calculated 
injections. Middle: One-minute sampling. Bottom: 30-minute 
sampling. 



Time from mid-transit (days) 

Fig. 15.— Delta (line) between a spherical and an oblate Jupiter¬ 
sized planet, using Saturn-like (/ = 0.1) oblateness, spin obliquity 
of 45° and a 87-day (Mercury) orbit. Data points are for 6 years 
of quiet solar observations using PLATO 2.0. No noise modeling. 


hope to measure the oblateness of hundreds of planets, 
and create an insightful distribution statistic. 

5. DISCUSSION 

After examining these transits, we should summarize 
the learnings. It is clear that future photometry will 
allow for many more discoveries, and we should prepare 
to get most out of these data. 

5.1. Observation strategies for PLATO 2.0: Single 
transits 

While transit detections of medium-sized, short-orbit 
planets such as Earth and Venus will be feasible with fu¬ 
ture photometry, the stellar noise will put limits to planet 
sizes and configurations for Mercury and Mars analogues. 
The greatest challenge will be longer-orbit transits (cold 
gas giant analogues), for which only one transit can be 
observed in a reasonable (<10 years) time span. With 
their long periods of 84 years (Uranus) and 165 years 
(Neptune), the observation of a single transit for any 
given solar system analogue is unlikely {Ptr = 0.024% 
for Uranus, Ptr = 0.014% for Neptune). The occurrence 
rate of such planets is currently not well constrained, as 
only a few examples have been found through gravita¬ 
tional microlensing (Furusawa et al. 2013; Sumi et al. 

2010), but is believed to be >16% at 90% confidence 
(Gould 2006). PLATO 2.0 will observe 85,000 bright 
(V<11) stars (Rauer et al. 2014), and all of these can 
be expected to yield high (S/N> 500) detection poten¬ 
tial for cold Uranus and Neptune analogues. For the 
full sample of >lbn stars, one might expect at least half 
of them to give sufficient S/N (> 20) for a clear detec¬ 
tion of such planets. Consequently, it can be expected 
to detect a considerable number of single-transits events 
for such planets with PLATO 2.0: With an occurrence 
rate of 16%, and 0.024% transit probability, and observ¬ 
ing 6 out of 84 orbital years, we can expect to find one 
Uranus analogue among 365,217 stars. For a 500mn sam¬ 
ple of PLATO 2.0 data, we can thus expect to find 1,369 
Uranuses (and not a single one among Kepler’s 150,000 
stars). 

The frequency of Jupiter analogues at distances of 
3-6 AU has been estimated as 5 ± 2% by Lineweaver & 
Grether (2003), or 3.3±1.4% (Wittenmeyer et al. 2011). 
Assuming a 5% frequency, we can expect to find ~5,000 
Jupiters and ~ 4, 000 Saturns, accounting for their higher 
transit probability (compared to Neptunes or Uranuses). 

It is clear that these numbers are only very rough es¬ 
timates, and might be subject to change by an order 
of magnitude. Through the very high number of stars 
observed with PLATO 2.0, however, it will be possible 
to observe single-transit events of all solar system ana¬ 
logues. Current exoplanet science relies on multiple tran¬ 
sits to determine periods, and confirm transit signals. It 
will be the challenge of these large-sample future mis¬ 
sions to treat single events in a way that allow for their 
detection and confirmation. This will put our own solar 
system in perspective: Is a configuration of rocky inner, 
and gaseous outer planets common, or exotic? 

5.2. Validation of transit signals 

There are three principle ways to validate transit-like 
signals. The first is to check if they are, at all, astrophys- 
ically possible. If this is the case, we can try to eliminate 
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those that are more likely to be caused by instrumental 
or stellar artifacts. 

The first check should be whether a putative star- 
planet-moon (or ring) configuration is astrophysically 
possible, plausible and stable over the long time. While 
the possibility of any configuration is essential, the situa¬ 
tion is less clear for the plausibility: Before the discovery 
of Hot Jupiters, such planets would have seemed implau¬ 
sible, and their discovery (Mayor & Queloz 1995) was 
challenged as they were found to be incompatible with 
theories of planetary formation (Rasio et al. 1996). 

This is also true for the stability. For example, we 
should assume that most moons are stable over long 
(Gyr) times, but there might be configurations for which 
this is not the case. A related example are Saturn’s rings, 
which are known to be unstable on timescales of <100 
Myr (Dougherty et al. 2009). Thus, instability should 
raise doubts about the presence of a moon (and less so 
about the presence of a ring), but does not proof their 
non-existence. 

To check the plausibility of moons, we can require that 
prograde moons are within ^39% of their Hill radius, and 
retrograde moons to be within 93% to be stable (Domin¬ 
gos et al. 2006). Also, they cannot have their orbit inside 
the Roche lobe, otherwise they would be torn apart and 
create a ring system. Similar arguments can be made for 
ring systems: When a flux increase (through diffractive 
forward-scattering) favours a ring composition of ice over 
rock, then such a ring planet cannot be too close to its 
host star, otherwise the ice would evaporate. Further¬ 
more, we can count the fraction of transits that exhibit 
no moon-like signals (assuming a coplanar moon), and 
compare this to the theoretical upper limit of ~6.4%, 
as derived by Heller (2012) using the Roche stability 
criterion. 

Sometimes, however, the situation is less clear. Re¬ 
cently, a single moon-like transit signal was reported by 
Cabrera et al. (2014) in the photometry of Kepler- 
90g. The best-fit parameters were a physically possible 
7.96 ± O.65i?0 planet orbited by a 1.88 ± O. 2 I.R 0 exo¬ 
moon. As reported by Kipping et al. (2015), the signal 
is well explained by an instrumental sensitivity drop of a 
single pixel on the CCD. The authors introduce a method 
dubbed “centroid map”, which compares the sensitivity 
of neighboring CCD pixels over time. This method al¬ 
lows to attribute probabilities to instrumental changes in 
sensitivities, for example caused by cosmic ray hits. 

While such sensitivity drops are severe, but usually 
short events, there might also be instrumental trends 
on timescales longer than the actual transit (s). It has 
proven effective to remove these by cross-correlating 
them with their physical originator (when the data is 
available), for example temperature drifts, or shifts on 
the CCD have been corrected with proxy data for Hubble 
Space Telescope data (Demory et al. 2015), and in the 
Kepler K2 mission (Vanderburg et al. 2015; Foreman- 
Mackey et al. 2015). 

When it comes to stellar noise, such proxy data is usu¬ 
ally not available, but could be obtained from observa¬ 
tions in multiple wavelengths through filters, or spectra. 
Stellar rotation can be measured (e.g. McQuillan et al. 
2015) and used to apply a suitable noise model. A few 
cases have been reported where transit depth variations 
might be explained by the stellar rotation phase which 


TABLE 1 

S/N FOR SOLAR SYSTEM TRANSIT OBJECTS, USING PLATO 2.0 
AND 6YRS OE quiet SUN DATA 


Object 

^ transits 

S/N 

Comment 

^ Mercury 

25 

7.0 


$ Venus 

9 

23.8 


© Earth 

6 

28.6 


([ Luna 

6 

2.1 


cf Mars 

3 

7.2 


\ Jupiter 

1 

4,360 


' ' Galilean moons 

1 

8.6 

out-of-transit 

b Saturn 

1 

3,575 


0 Saturn’s rings 

1 

21.5 

out-of-transit 

5 Uranus 

1 

714 


2 Neptune 

1 

783 



causes the transiting planet to occult star spots on a 
highly spotted star (Croll et al. 2015). With a com¬ 
bined detection of TTVs and rotation, it is also possible 
to distinguish between prograde and retrograde motion 
(Holczer et al. 2015). 

These measurements make actual use of the otherwise 
unwanted stellar variation. In order to only distinguish 
between stellar flux variation and a transit signal, sta¬ 
tistical tools can be used. We suggest to compare the 
occurrence rate of the signal(s) in question, including 
their shape, to the whole dataset in suitably chosen time 
bins. Ideally, there should be no other features in the 
whole dataset compared to the one in question. To vali¬ 
date exomoon signals, Hippke (2015) argued that a test 
should be made, checking whether it is possible to shift 
the (virtual) folded transit time (to any other position), 
and then still have a significant dip. A variation of this 
test can check the uniqueness of any given signal. 

5.3. Data handling burdens 

Given the required detection techniques, it will be an¬ 
other challenge to efficiently mine the large data volume 
expected by PLATO 2.0. When extrapolating the Ke¬ 
pler data volume (^16 MByte per star and year in Imin 
integrations), we will have to search 48 PByte (48 x 10^^ 
Byte) of data. The storage of these data alone costs 
^480,000 USD, in 2015 storage prices. Hard disk capac¬ 
ities have increased by a factor of 16 (for constant nomi¬ 
nal prices)^ between 2004 and 2014. Extrapolating to the 
PLATO 2.0 mission in 2024 will estimate storage require¬ 
ments of ^30,000 USD. Clearly, most analyses will have 
to be done remotely, with the data stored in central facil¬ 
ities such as the Mikulski Archive for Space Telescopes 
(MAST)®. Searches in Kepler data today can be done 
on the researcher’s own machine, as the total volume of 
the dataset is < 1TB, i.e. less than 100 USD in 2015 
prices. Depending on the computational demands of a 
specific research, processing is done on single computers. 
University Clusters, or the NASA Supercomputer. We 
expect that most or all large-scale searches in PLATO 
2.0 data will need to be performed remotely, using super¬ 
computer facilities. 

5.4. Planets interior to Mercury 

Our solar system does not have any planets interior to 
Mercury, but this is not the norm: It is estimated that 

® http://www.jcinit.com/diskprice.htni, retrieved 09-Mar 2015 

® https://archive.stsci.edu/ 
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Fig. 16.— Period-radius parameter space for planet transit de¬ 
tectability. The threshold (line) is discussed in section 5.4. 

^5% of stars host (multiple) tightly-packed inner planets 
(Lissauer et al. 2011), and ~50% of stars have at least 
one 0.8 — 2i?0 planet within Mercury’s orbit (Fressin et 
al. 2013). The latest simulations indicate that our solar 
system also possessed such planets, but these have been 
destroyed in catastrophic collisions, leaving only Mer¬ 
cury (Volk & Gladman 2015). To better characterize 
these configurations of solar system analogues, it would 
be helpful to be able to detect them in future missions. 
As we will see in the next section, we can expect to find 
most or all planets > O.5i?0 within Mercury-sized or¬ 
bits with PLATO 2.0, answering interesting questions 
on planetary formation. 

5.5. Ultimate limits of photometry 

To determine the actual limits of detection in the ra- 
dius/period parameter space, we have created a diagram 
of achievable S/N using a series of injections and blind 
retrievals. We used our Sun’s data with its red noise, 
and injected series of transits on circular orbits with 
central transits, i.e. impact parameter 5 = 0. Limb 
darkening parameters were hxed (“known”) for simplic¬ 
ity, as we found that their typical uncertainties account 
for < 1% of transit depth/duration variation. To check 
this, we have used typical uncertainties in stellar parame¬ 
ters and recalculated the transits with varying limb dark¬ 
ening coefficients from (Claret & Bloemen 2011). For 
example, in the recent characterization of Kepler-138 by 
Jontof-Hutter et al. (2015), the stellar parameters were 
given as T^ff = 3841 ± 49K, [Fe/H]=0.280±0.099, and 
log(g)=4.886±0.055. With these uncertainties, the tran¬ 
sit depth varies by < 0.5%. 

We only varied the radius of the injected planet, and 
the period. Retrieval was done blindly, and considered a 
success if both parameters were within 10% of the true 
values. 

We required the likelihood of a noise-based dip of the 
same depth and shape as the transits to be < 1 : 1,000 in 
6yrs of data, using blind recovery and our noise model, 
plus the PLATO 2.0 instrumental noise. With this limit, 
one false positive among 1,000 real planet detections is 
found. The measured S/N as described in section 2.4 for 
such a borderline detection is ~10. If the noise model 
is neglected, one false positive among 300 real detec¬ 
tions occurs; to achieve the same false positive rate, the 


threshold would need to be increased to S/N~14. The 
result in Figure 16 shows that planets > O.5i?0 within 
Mercury-sized orbits can be detected. What is more, we 
can expect to detect all planets in the habitable zone of 
a quiet G2-dwarf, i.e. on 200-600 days periods and radii 
> O.7i?0. This result can easily be transferred to other 
stellar radii. Given sufficient brightness and the same 
low stellar noise, the transit depth scales as /i?*. For 
example, a Mars-sized planet of O.53i?0 around a 0.5Rq 
M-dwarf would cause a transit dip of ^97ppm, which 
is about the depth of Earth’s transit orbiting our Sun. 
Consequently, and helped by shorter periods (more tran¬ 
sits), we can expect to detect all potentially habitable 
planets orbiting M, K and G stars with 6yrs of PLATO 
2.0 data. 

It is also worth noting that the detectability function 
in Figure 16 has a step at P > l/2Tobs (not shown in 
the Figure), i.e. when only one transit is observed during 
the mission duration. Then, all else equal, longer periods 
are preferred due to their longer transit duration which 
makes the detection easier. For comparison, Neptune 
has a transit duration of 3.Id, and Jupiter “only” 1.4d, 
giving a S/N advantage of a/3.1/1.4 ^ 1.5x from transit 
duration. Of course, the occurrence rate of Neptunes and 
Jupiters is likely different (see section 5.1), as is their 
transit probability. It is important, however, that the 
S/N differences are taken into account when calculating 
population statistics with PLATO 2.0 data. 

On the other side, to answer the question of what near¬ 
perfect photometry can deliver in the future, we have 
repeated our injections without the instrumental noise. 
The results for the G2 star are very similar, as most of the 
noise already comes from stellar trends in the PLATO 2.0 
scenario. As has been shown in section 4.2, the improve¬ 
ments towards PERFECT are only ~10%, and up to 
^50% for the smaller M-dwarfs. According to our blind 
retrieval simulations, even with the inclusion of standard 
noise models, a detection of Mercury or Mars analogues 
is beyond the photometric limits. This, however, is not 
necessarily the ultimate limit of photometry (as asked in 
the introduction); we might improve noise models using 
proxy data from radial velocity, spectroscopy, hydrody- 
namical modeling, and other methods yet to be devel¬ 
oped. 

6. CONCLUSION 

In this work, we have shown that future photome¬ 
try will be able to detect Earth- and Venus-analogues 
when transiting G-dwarfs like our Sun (Table 1). Larger 
sized planets (> 2 P 0 ) will be detected in a single transit 
around G-dwarfs, in low stellar noise cases, and assum¬ 
ing one can find them in the Hrst place. The search 
techniques for such single transits will require further 
research and validation, and will likely be performed re¬ 
motely, due to the large storage requirements. 

While the detection of moons in a solar system config¬ 
uration will remain problematic in the next decades, the 
situation is better for rings, and for moons in M-dwarf 
systems. 

For source stars with strong red noise characteristics, 
such as our Sun, we suggest to shift the usual S/N limit 
from 7 to 14, in order to prevent too many false posi¬ 
tives. This limit can be set individually for each star in 
question, by multiply injecting and retrieving artificial 
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signals to the one that is in question. Alternatively, we 
recommend to model the noise and set an appropriate 
threshold for the noise model using simulations. We re¬ 
lease all data used in this paper, plain and injected, for 
the community' and encourage testing different (and also 
blind) retrieval techniques in preparation for PLATO 2.0 
and other missions. 

Despite these challenges, we believe that the era of 
transit planet detection is still in its infancy, and pho¬ 


tometry will have a bright future in the coming decades. 
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